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NOMENCLATURE 






r/ 


A = flow area, 

d = diameter, m 

F = thrust, N 

g = acceleration due to gravity, 9.8 m/s^ 

M = Mach number 

P = static pressure. Pa 

Pj = stagnation pressure. Pa 

Q = heat addition rate, w 

q = heat addition per unit mass, J/Kg 

T = static temperature, °K 

Tj = stagnation temperature, °K 

V = velocity, m/s 

6 = conical diffuser half-angle, degrees 

Y = ratio of specific heats 

(j = density, Kg/m^ 

o = acceleration, g's 

Ki, A, B, C, D = constants used in Eqs. 8, 9, 10, and 11 

A numerical subscript on a variable corresponds to a fluid state 
defined in Figure 3a and/or 3b. 


iv 


INTRODUCTION 


Direct launch systems are currently being considered for space re- 
supply and disposal missions where high accelerations can be tolerated. 
Payloads delivered into space with these systems would be accelerated 
to very high velocities (on the order of 10'^ m/sec) close to the earth's 
surface and they would then be slowed somewhat to the desired final 
energy during passage through the earth's atmosphere. One such system 
that has been proposed recently for such missions is the electrothermal 
ramjet.^ This system appears to be capable of achieving the high ac- 
celeration levels needed to assure launch tracks that are of a reason- 
able length primarily because it does not accelerate either the 
propellant or the power plant required to effect its operation. 

Conceptually the electrothermal ramjet involves a long tube filled 
with gaseous propellant through which a ramjet engine is accelerated to 
a specified terminal velocity. Typically the ramjet would be operated 
at constant thrust using hydrogen propellant heated by either electro- 
magnetic radiation or electrical currents as it passes through the 
ramjet. Energy would either be beamed or conducted to the ramjet heat 
addition zone from the tube walls as the ramjet passed down the launch 
tube. Because a ramjet requires a substantial initial velocity before 
it can begin to produce thrust it is anticipated that the ramjet would 
be accelerated to this velocity by a device such as a light gas gun. 
Preliminary analysis ’ has suggested that a particular ramjet design 
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referred to as the annular flow, plug ramjet is well suited to this 
application and that such a device should exhibit good performance. 

The general concept of the electrothermal plug ramjet is shown 
in Fig. 1. The projectile, which is a solid centerbody, is shown 
traveling down the circular tube. The reference frame is chosen to be 
moving with the projectile so the gaseous propellant is shown moving 
towards the projectile with velocity V. The incoming propellant will 
encounter the annular diffuser where it will be slowed relative to the 
projectile. The coupling of the tube wall with the projectile forebody 
will produce the geometric configuration necessary to bring the propel- 
lant to subsonic conditions at the exit of the annular diffuser zone. 
After the propellant has been brought to subsonic conditions, heat is 
added in the annular heat addition zone, and the propellant then passes 
through the annular nozzle zone, which is a converging-diverging chan- 
nel, and is expelled out the rear of the device so as to produce a 
thrust force on the projectile. This annular design is utilized be- 
cause it has several attractive features including 

1) Heat addition can be effected outside the ramjet body directly 
from the tube wall. This eliminates the need for windows in 
the ramjet body for the case of beamed energy heating and 
sliding current contacts for arc discharge heating schemes. 

2) Direct mechanical or fluid mechanical forces induced between 
the tube wall and the projectile should tend to stabilize 
and center it. 

3) The structural design of the projectile appears to be rela- 
tively simple. 

4) Shock losses in tfie diffuser can be minimized by properly 


designing the annular diffuser section. 

5) There will be no "spill over" of propellant. 

6) Variations in fluid flow areas needed to achieve proper flow 

conditions as the projectile velocity increases can be 

achieved by varying the cross-sectional area of the launch 

tube as a function of axial position. 

Preliminary analyses regarding annular flow electrothermal ramjet 

1 2 

operation have been conducted ’ and typical results obtained from these 
studies showing design and operating parameters for a typical mission 
are presented in Figure 2. These results assume that there are no 
diffuser shock losses, that the hydrogen propellant behaves as a perfect 
gas, and that the fluid flow is quasi-steady. These results detail the 
performance of a ramjet launch cycle for the following specific opera- 
tional parameters 

1) A constant projectile mass of 10 kg. 

2 ) A constant acceleration of 30,000 g's (hence constant thrust) 
from rest to 15,000 m/s. 

3) An initial propellant temperature of 300°K. 

4) An initial propellant pressure of 30 atm. 

5) A fixed projectile diffuser throat diameter of 16 cm (Fig. 1). 

6) A fixed projectile heat addition zone diameter of 8 cm (Fig. 1). 
If one assumes that the projectile will be acted upon by either the 

ramjet force or some externally applied force needed to produce the de- 
sired acceleration, the projectile trajectory will be known. With the 
projectile velocity profile known, the propellant flow velocity relative 
to the ramjet is known at each position along the launch track. One can 
now prescribe the diffuser throat area required at each position along 
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the tube so that the flow will remain sonic at this throat. One can 
also prescribe the heat supplied in the subsonic heat addition region 
so that the specified thrust will either be produced or the Mach number 
at the exit of the heating region (in the projectile frame of reference) 
reaches its maximum allowable value of unity. After heat has been added, 
the flow is accelerated isentropical ly in the nozzle to the maximum 
supersonic Mach number achievable under the constraint of the area 
change available. The static pressure at the nozzle exit is determined 
by this expansion and it will vary with projectile position. Proper 
control of this pressure at the base of the projectile would have to be 
achieved in an actual device by proper design of the launch tube. Such 
a tube would consist of two sections: a cavity section into which 

hydrogen processed by the ramjet would be collected, and the launch sec- 
tion through which the projectile would pass. It has been assumed in 
this analysis that the cavity section could be designed to have the 
length and cross-sectional area needed to accommodate the propellant and 
facilitate maintenance of the necessary projectile base pressure. While 
it seems possible to do this intuitively, the examination of the unsteady 
fluid mechanics describing the process has not yet been attempted. 

Figure 2 shows the results of such an analysis for the case where 
a constant thrust of 3 x 10^ N is maintained on a 10 kg mass. As Fig. 2a 
shows, the projectile velocity (in the launch tube frame of reference) 
increases linearly with time. Since the propellant is initially at rest 
relative to the tube wall the projectile velocity is equal to the pro- 
pellant inlet velocity in the projectile frame of reference. The ex- 
haust velocity (in the projectile frame of reference) is observed to be 
equal to the inlet velocity during the first -v 8 ms of operation, so 
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during this time interval, essentially no net thrust is produced through 
the mechanism of momentum change of the propellant. After about 16 nis 
of operation, however, the exhaust velocity has risen to about 15 km/s 
and has begun to level off. During the time between 8 and 16 ms, 
all of the heat that can be added without choking the flow {achieving 
unity Mach number at the end of the heat addition region) is being 
added. 

As Fig. 2b shows, the net ramjet thrust is inadequate to meet the 
3 X 10^ N thrust requirement at times less than 16 ms. Thus an alterna- 

3 

tive accelerating scheme, such as the light gas gun would have to be 
used during this time interval. As the projectile velocity increases, 
more heat can be added and this causes the increase in exhaust velocity 
shown in Fig. 2a. After 16 ms, both the exhaust velocity and the mass 
flow rate have risen to the point where the desired thrust level can be 
produced and the Mach number after heat addition begins to drop below 
unity. 

Figure 2c, which indicates projectile position as a function of 
time for this case of constant acceleration, shows a launch tube about 
400 m long would be required to accelerate the 10 kg mass to the ter- 
minal velocity of 15,000 ms in 'v. 50 ms. Figure 2d shows the variation 
in tube and projectile diameters that would have to be effected to sat- 
isfy fluid mechanical constraints imposed on the ramjet. As previously 
stipulated, the projectile diameters at the diffuser throat and in the heat 
addition region have been treated as fixed at 16 and 8 cm respectively. 
The reduction in tube diameter required to yield unity Mach number at 
the diffuser throat as the projectile velocity increases is shown by 
the solid line. Data from this line taken together with the position 
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profile of Fig. 2c specify how the tube diameter would have to vary 
along its length. Figure 2d also shows that the projectile diameter at 
the nozzle throat would have to increase as the projectile velocity in- 
creases and the input power is varied to induce the desired thrust level. 
The fact that this diameter must change during the launch cycle repre- 
sents a significant mechanical complication in the system. It is noted, 
however, that the amount of change In this diameter becomes very small 
after about 20 ms when the projectile velocity is 'v 5 km/s. It is an- 
ticipated that a pressure- or temperature-sensitive scheme might be 
designed to facilitate this change in projectile diameter at the nozzle 
throat. Alternatively, one could hold the projectile diameter in the 
nozzle region fixed, reduce the tube diameter slightly more rapidly than 
Fig. 2d suggests, and abrade the diameter of the projectile at the dif- 
fuser throat using protrusions from the tube wall. 

Figures 2e, 2f, and 2g indicate how the stagnation temperature in 
the heat addition zone and the static pressure and temperature at the 
nozzle exit should vary to conform to the problem specifications. The 
peaks in temperature and pressure occur at the point where the power 
input is a maximum (Fig. 2h). This, in turn, occurs at the time when 
the thrust has just reached its design value and the Mach number at the 
end of the heat addition region is still at unity. At this condition, 
the exhaust velocity is still relatively low and the bulk of the thrust 
is being produced by the static pressure difference across the projec- 
tile. If one did not start using the ramjet until after 'v 20 ms into 
the launch, the exhaust pressure and temperature would lie in a very 
moderate range throughout the period of ramjet operation. Note that it 
is the pressure profile of Fig. 2f that would have to be matched 
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through proper design of the cavity section of the launch tube. This 
might be achieved by selecting not only the length and cross-sectional 
area of this section, but also by the inclusion of accumulator chambers 
that could be connected into the tube at the proper time. 

The heat addition zone stagnation temperature (Fig. 2e) is ob- 
served to remain near 10,000°K after about 16 ms. It is noted that 
dissociation and ionization may introduce significant effects while 
operating within this regime, and these effects will be examined later 
in this analysis. 

Figure 2h shows the thermal power input profile required to effect 
operation as well as the actual profile of mechanical power input to the 
projectile. After the 20 ms period over which the auxiliary starting 
force would be applied, it is observed that the ratio of projectile 
power to thermal input power rises above -v- 50%. The ramjet is observed 
to become increasingly efficient as the launch cycle proceeds. This 
suggests that the ramjet will be more efficient at the higher velocities, 
where the efficiencies of electromagnetic launchers drop off. The dif- 
ference between the projectile and thermal input powers represents the 
power being deposited as kinetic and thermal energy in this hydrogen 
propellant. It should be noted here that frozen flow losses in the 
nozzle and electrical-to-thermal power conversion losses as well as 
frictional losses have been neglected in this preliminary analysis. 

The results presented in Fig. 2 suggest that the annular flow 
electrothermal ramjet concept has the potential of achieving efficient 
high velocity high acceleration launches, and also that this concept 
cannot be rejected on the basis that it violates a basic principle of 
fluid mechanics or thermodynamics. The pressure, temperature and power 
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demand profiles are all reasonable, and it is for these reasons that 
further analysis of the annular flow electrothermal ramjet was under- 
taken to bring to light aspects which may degrade ramjet performance. 

There are several effects which may degrade plug ramjet performance 
below that shown in Figure 2. Some of these effects are 

1) Diffuser shock losses. 

2) Propellant dissociation, ionization and non-perfect modes of 
molecular and atomic energy storage. 

3) Freezing (translational, rotational, vibrational , chemical). 

4) Non-steady gas dynamics effects. 

It is the purpose of this thesis to address these effects and to de- 
termine the magnitude of degradation of annular flow electrothermal 
ramjet performance that they might cause. 

In order to illustrate the performance of the electrothermal plug 
ramjet, an analysis will be also conducted on the more familiar con- 
ventional ramjet which uses hydrogen propellant that behaves as a 
perfect gas. The results from this analysis will be compared to analyses 
conducted for the annular flow electrothermal plug ramjet. Simple 
theoretical models detailing conventional ramjet and electrothermal 
ramjet operation with perfect gas hydrogen propellant will be presented, 
and data from the more extensive analyses considering dissociation, 
ionization and non-equilibrium processes will be interpreted in the 
light of the results from the simplified models. 

It should be stated here it has been assumed in all analyses that 
the fluid flow is axis -symmetric or one-dimensional and quasi-steady, 
and the fact that the tube wall is moving in the reference frame of 
the annular flow electrotheimal plug ramjet is assumed to have no 


influence on the flow. All flowfields analyzed are considered to 
be inviscid. The ideal gas equation of state is assumed to be valid 
at all times. 



FLOWFIELD MODELS 

The general concept of a conventional ramjet is shown in Fig. 3a. 
The reference frame is chosen to be moving with the ramjet, so the pro- 
pellant is shown moving towards the inlet with a velocity . The Mach 
number associated with the incoming flow is designated M^ . After the 
ramjet has been accelerated to a high enough initial velocity, an ob- 
lique shock will be attached to the conical tip of the center body 
which serves as a diffuser. Downstream of the diffuser tip, a series 
of shocks reflected between the center and outer bodies will occur. 

They will terminate in a normal shock that will slow the flow to sub- 

4 

sonic conditions. A review of the hypersonic wind tunnel literature 
has suggested that normal shock recovery can be expected after the 
initial conical shock, so it is assumed that the flow passes through an 
initial conical shock to State 2 and then passes through a single nor- 
mal shock to state 3. One-dimensional values obtained by averaging 
over the axis-symmetric conical flow field variables computed at the 
exit to the conical flow field in Region 2 are used as input for the 
normal shock relations. For a description of the methods used to com- 
pute the average conical flowfield values see Appendix A. 

Across the normal shock, equations describing mass continuity, the 
first law of thermodynamics, conservation of momentum, and the second 
law of theniiodynamics are satisfied. Immediately downstream of the 
normal shock it is assured that the area between the center and outer 
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bodies is sufficient so the flow will pass through the minimum area 
that is encountered at the diffuser throat. After the normal shock, 
the flow enters the subsonic heat addition region at 4 and heat is 
assumed to be deposited into the flow between States 4 and 5. In the 
heat addition region, the first and second laws of thermodynamics are 
satisfied, along with mass continuity and the momentum equation. After 
leaving the heat addition region at 5, the flow passes through a 
converging-diverging channel where it is accelerated first to the sonic 
condition at state 6 and then on to supersonic conditions at the nozzle 
exit 7. /c? 



With this conventional ramjet^^igure 3a) the static inlet and 
outlet pressures are equal becaus^the free stream conditions are com- 
municated outside of the moving outer body. In order to satisfy this 
condition while assuming no shocks in the nozzle one finds the nozzle 
exit area (Ay) must be larger than the free-flow diffuser inlet area 
(A-j). It is hoped that by comparing this more familiar ramjet to the 
new one being proposed that the concept and performance of the plug 
ramjet can be better understood. 

The annular flow plug ramjet is shown in Fig. 3b. It appears 
similar to the conventional ramjet at first glance, but is different in 
the following ways 

1) The plug ramjet moves down the tube as the center body of the 
conventional ramjet does but in the plug ramjet case the outer 
body is the launch tube which is stationary. Because the 
launch tube diameter is essentially constant along its length 
the fluid mechanics and thermodynamics of the problem can be 
considered from the frame of reference of the ramjet as though 


this tube were not moving. 

2) The inlet and outlet areas associated with flow through the 
plug ramjet (A^ and A^) must be equal, again because the launch 
tube diameter does not vary significantly along its length. 

3) The inlet and outlet pressures and are not equal be- 
cause the upstream and downstream regions are not connected 
through an exterior path. In fact the outlet pressure P^ will 
be greater than the inlet pressure P^ and this pressure dif- 
ference will therefore produce some component of the thrust 
exerted on the plug. 

One can see that the fluid mechanics of both configurations in 
Figs. 3a and 3b are described by the same basic equations. The only 
difference in the two cases is the difference in the boundary condi- 
tions that are applied to the flow. 

When non-perfect gas modes of energy storage such as ionization 
and dissociation are allowed to occur, results from statistical mechan- 
ics are utilized to compute equilibrium compositions and thermodynamic 
properties and shifting chemical composition is allowed to occur. For 
a description of the methods used to compute the equilibrium properties 
of hydrogen see Appendix B. Where they are deemed potentially impor- 
tant, nonequilibrium processes are evaluated by considering finite rate 
reactions and by calculating the conditions at which specific reactions 
will become slow enough so the chemical composition of the flow becomes 
essentially frozen, A description of the methods utilized to perform 
these calculations can be found in Appendix C. 

Because the essense of the phenomena occurring in ramjet devices 
of the types considered here can be described by simple, perfect gas 
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models, these types of theoretical models will be presented to illus- 
trate the basic operational trends present in the analyses conducted. 
Perfect Gas Models 

Consider now that either ramjet in Fig. 3 is being accelerated at 
a constant thrust F and that heat is being added at a rate Q to achieve 
this level of thrust. In general Q will vary with the velocity of the 
ramjet and will be different for the conventional and plug designs. By 
assuming that the propellant flowing through the control volumes shown 
in Figs. 3a and 3b behaves as a perfect gas, the momentum equation ap- 
plied across the control volumes of Fig. 3 can be put in the following 


T, 1+ M? 


Y-1 M2l Y-1 


p-f - = • 


2 "7 


In this equation the entrance and exit conditions are respectively; 
Mach number M-j , M^; area k-j and stagnation pressure , Pj^. The 
ratio of specific heats is y and P^ represents the entrance static 
pressure. The left hand side of this equation represents a non-dimen- 
sional ized thrust. The stagnation pressure ratio appearing in Eq. 1 


can be written 


i A J5 _I7 

D ” P p p 

^4 ^5 


where (Py /P^. ), the stagnation pressure ratio across the diffuser, is 
'4 'l 

less than unity because of diffuser losses. The heat addition region 

stagnation pressure ratio (P, /?^ ) is also less than unity because of 

'5 '4 

heat addition-induced stagnation pressure losses. Because the nozzle 

is assumed to be isentropic the ratio (9^ /P^ ) is unity. 

‘7 ‘5 





When the first law of thermodynamics is applied to a perfect gas 
working fluid, the following expression for the instantaneous ramjet 
efficiency (power into ramjet mass acceleration divided by input thermal 
power) is obtained 



F Yi 7.,, 


(3) 



/ 




To apply these equations to the analysis of either ramjet it is neces- 
sary to specify the inlet Mach number (M-j) and the diffuser geometry so 
that the diffuser stagnation pressure ratio (?, /Py ) can be computed. 

‘4 '1 

For a given conical diffuser half-angle (6) and annular diffuser shoul- 
der diameter (d^) and heat addition center body diameter (d^) the inlet 
free-flow diameter (d^) is specified by placing the intersection of the 
conical shock and the outer body or tube wall at the annular diffuser 
shoulder. Specifying these parameters serves to fix the flow conditions 
up to State 4, where subsonic conditions are always specified. 

The flow downstream of the conical flowfield is assumed to be one- 
dimensional, so one-dimensional values for the conical flowfield var- 
iables are computed by the methods detailed in Appendix A, These one- 
dimensional values of the variables are then used as input to the nor- 
mal shock relations so the stagnation pressure ratio /P-r and the 

I4 , 

flow conditions downstream of the normal shock at 3 can be utilized to 
specify the conditions at the entrance to the heat addition region 4. 

By knowing and specifying M^, the heat addition stagnation pressure 
ratio (Py /P-r ) is fixed, as is the stagnation temperature ratio 

'5 ‘4 

(Tt /Ty = T-. /T, ). To analyze ramjet operation then, either the 

'5 '4 '7 '1 
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exit/ inlet area ratio (A^/A^) or the exit/inlet static pressure ratio 
(Py/P-l) may be specified to fix For the conventional ramjet the 
pressure ratio Py/P^ is set to unity, and for the electrothermal ramjet 
the area ratio A^/A^ is set to unity. It is noteworthy that the param- 
eters A^/Ap P^/P^ and are not completely independent. Specifying 
one parameter dictates values for the others. 

Dictating the exit/inlet conditions for the ramjet control volume 
allows the calculation of other parameters such as the exit/inlet velo- 
city and temperature ratios, and the fractions of the total thrust pro- 
vided by 1) the difference in the pressure area product between the 
entrance and exit flowfields and 2) the momentum change in the propel- 
lant. The mathematical relationships utilized to compute these ratios 
and fractions are not shown, but graphical results are illustrated in 
the Results Section of this thesis for specific operating conditions. 

To further illustrate the performance of the ramjets operating 
with a perfect gas propellant it is instructive to express the first 
law of thermodynamics in the form 


1 \\/ 


/ ^7 Pi 


and the momentum equation in the form 


!z 1^ + 1 

h " 7 "i^ ^T^ Vt^ 7 


TvTvfT 


Equation 5 can be multiplied by the quantity 


0 



_J_ D T 

Y-1 


to obtain the form 


('^Ofe i '' 


1 \ '1 


For the conventional ramjet 


P7/P1 = 1 . 




The stagnation pressure ratio P^; /P_ is given by 

'7 '1 


’’t ^ ^ 

-_i = Ji 'y 

Pf Py Py Py 

'1 '1 '4 ^5 

where P, /P, is fixed by the inlet flow conditions and the diffuser 
4 1 

geometry and P^ /P^ is unity. The stagnation pressure ratio P,. /P^ 
'..75 '5 '4 

is dictated by the rate of heat addition per unit mass to the propellant 

and for subsonic heat addition this ratio does not vary appreciably. 

Examination of Eq. 4 shows that for fixed values of M, , y, /P^ , P./P, 

T4 Ti 1 7 

and Pj /Pj , the effect of the variation of P^ /P^ on the expression 
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is small compared to the change in 
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as the rate of heat addition per unit mass is varied. As a result, 
Eq. 4 can be written 



where is approximately constant. The variation in is typically 
of the order of 5% for the cases considered. 


Optimum Ramjet Operation 

It will be shown in the Results Section of this thesis that there 
is a value of the non-dimensional thrust [F/P^A^] for which the instan- 
taneous thermal efficiency [FV^/Q] is at a maximum for a specified 
diffuser and given free-flow inlet conditions. An approximate mathe- 
matical relationship which predicts the point of maximum efficiency for 
the conventional ramjet operating with a perfect gas propellant can be 
derived by considering the first law of thermodynamics and the momentum 
equation. The quantity appearing in Eq. 5 



varies by about 6.5% at the most, so it can be treated as a constant. 

By assuming that the previously mentioned quantities are constant Eq. 5 
can be put into the following functional form 
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where A and B are constants. This equation states that the ramjet 
thrust varies as the square root of the heat addition rate per unit 
mass. Similarly the expression for instantaneous efficiency (Eq. 7) 
can be put into the form 



where C and D are constants. 

Examination of Eq. 10 reveals that there is an optimum stagnation 
temperature ratio that will yield a maximum instantaneous efficiency. 

From Eq. 9 one sees that this optimum temperature ratio corresponds to 
an optimum normalized thrust level. This optimum corresponds to a ram- 
jet characterized by specified geometry and inlet conditions. Changes 
in the inlet condition will necessitate changes in the optimum stagnation 
temperature ratio and the optimum normalized thrust level. The values 
for the constants C and 0 are given by 


C = K 2 (K^ + K 3 K^) 
D = K^Cl + K 3 ) 


where 
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These relations are valid where P^/P^ is held to be constant and the 

stagnation pressure ratio P- /P^. takes on typical values in the range 

'5 '4 

0.8 to 1.0. By differentiating the instantaneous efficiency (Eq. 10) 
with respect to the stagnation temperature ratio, the stagnation temper- 
ature ratio which gives the optimum instantaneous efficiency can be 
computed. The expression for this optimum ratio is 



Q 

Dissociating Flow Mode l 

It was anticipated that assuming the propellant behaves as a per- 
fect gas would introduce errors into these analyses, because it was 
recognized that effects such as vibrational and electronic excitation dis- 
sociation and ionization of the propellant would play significant roles 
in the performance of the ramjet device. The annular flow electrothermal 
plug ramjet was also analyzed with the above effects included. The va- 
lidity of assuming either chemically frozen or equilibrium chemical com- 
position in the supersonic nozzle expansion was also evaluated. In 
analyzing the flow in the diffuser for this dissociating flow case it 
was assumed that the fluid in the conical flowfield between the normal 
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shock and the oblique shock behaved as a perfect gas and thus that this 
portion of the model developed for the perfect gas models could be 

applied. This assumption should not introduce too much error into this 

analysis because the oblique shock angle should be shallow for high in- 
let velocities and small diffuser cone angles, and the temperature rise 

in the conical flowfield should not be too large. Downstream of the nor- 

mal shock, chemical relaxation and other non-perfect gas modes of energy 
storage were allowed to occur in the propellant. The propellant was as- 
sumed to be in thermodynamic equilibrium downstream of the normal shock and 
up to the annular nozzle throat (6). Finite chemical reaction rates were 
evaluated in the supersonic nozzle expansion between 6 and 7. Because 
shifting chemical equil ibrium was allowed, closed form equations for 
the thermodynamic properties could not be obtained. The method of 
analysis utilized for the model including shifting composition involved 
piece-wise application of the momentum and energy equations between 
each state identified in Fig. 3b that is downstream of the normal shock. 

In order to compare the performance of the perfect gas ramjet to 
that of the ramjet utilizing the chemically reacting propellant, the 
diffuser operating conditions are specified and the non-dimensional 
thrust is computed along with various other non-dimensional parameters 
such as the inlet/exit temperature and velocity ratios. Different 
values for the non-dimensional thrust [F/P^A^] are determined by in- 
dependently varying each component in the non-dimensional thrust equa- 
tion, and the downstream flow conditions are found such that a speci- 
fied, constant thrust level is met. 



RESULTS 


It is Instructive to examine a ramjet operating at a typical 
flight velocity that would be encountered at one point in a launch se- 
quence so that the influences of the various operational parameters 
can be established. In order to do this several different ramjet con- 
figurations have been analyzed for a moderate flight velocity of 
7000 m/s. A diffuser cone half-angle (6) of 30°, a static temperature 
(T^) of 300°K, and diffuser shoulder diameter-to-tube diameter (d^/d^) 
and heat addition region diameter-to-tube diameter (d^/di) ratios of 
0.805 and 0.403 respectively (Fig. 3), have been used in the analysis. 

The 30° cone angle assured conical shock attachment and reasonable diffuser 
efficiency. The static temperature T^ was selected as a reasonable 
value that could exist in an actual device. The diameter ratios re- 
sulted in adequate center body structural integrity and they were also 
such that choking at the diffuser throat was precluded. The specific 
ramjet configurations examined are (1) the conventional ramjet utilizing 
diatomic hydrogen which behaves as a perfect gas as a working fluid and 
with the constraint that the inlet and outlet static pressures are equal, 
(2) the annular flow electrothermal ramjet operating with the same per- 
fect gas and with the constraint that the inlet and outlet flow areas 
are equal and (3) the annular flow electrothermal ramjet operating with 
hydrogen propellant that is allowed to dissociate and ionize and to have 
vibrational and electronic modes of energy storage along with finite 



chemical relaxation rates occurring in the supersonic nozzle expansion. 
For this latter case two thrust levels per unit cross sectional area of 
tube have been specified. The higher value for the thrust per unit area 
is F/A^ = 9.5 X 10'^ Pa and the lower value for the thrust per unit area 
is F/A^ = 9.5 X 10^ Pa. Two values have been considered so that the 
effect of varying the initial pressure in the launch tube can be illus- 
trated for high and low thrust per unit area values. 

Figure 4 illustrates the relationship between the instantaneous 
efficiency [FV^/Q] and the non-dimensional thrust [F/P^A^] for the 
previously mentioned operational parameters. Increasing the non- 
dimensional thrust implies increasing the heat addition per unit mass to 
the propellant. The perfect gas conventional ramjet is shown to exhibit 
the best performance, and the optimum value for instantaneous efficiency 
is seen to occur at a non-dimensional thrust of 12 for the conventional 
ramjet. This peak value is predicted to first order by Eq. 10. 

The consequence of not expanding the supersonic nozzle flow to the 
pressure ratio Py/P^ of unity is illustrated in Fig. 4 by the curve for 
the electrothermal plug ramjet operating with a perfect gas (no dis- 
sociation). Its performance is seen to be degraded substantially with 
respect to the conventional ramjet operating with a perfect gas. The 
other curves, which illustrate the effects of chemical freezing, disso- 
ciation, ionization and vibrational and electronic modes of energy 
storage, are shown to exhibit performance which is similar but inferior 
to the electrothermal plug ramjet operating with a perfect gas. It is 
noted that the dissociating propellant curves of Fig. 4 were generated 
by fixing the thrust-to-cross sectional area ratio (F/A-|) and varying 
the initial pressure in the launch tube P-j . The maximum value of the 
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non-dimensional thrust for which each curve shown in Fig. 4 terminates 
is that value for which no more hea': could be added in the heat 
addition region of the respective ramjets without choking the flow at 
the exit of the heat addition region 5. All of the curves in Fig. 4 
illustrate, that for a given velocity and diffuser geometry, there is a 
value for the non-dimensional thrust which provides the maximum in- 
stantaneous efficiency. It is also evident that, for the case of a 
dissociating propellant, a high tube pressure high thrust device cor- 
responding to the higher thrusv. per unit area (F/A-j = 9.5 x 10^ Pa) 
provides performance superior to a low tube pressure low thrust device 
(F/A-j = 9.5 x 10^ Pa). This is due to the greater level of dissocia- 
tion that occurs at the lower pressures. Figure 5 illustrates this by 
showing that the mass fraction of diatomic hydrogen in the ramjet ex- 
haust is greater at the higher thrust per unit area level '■■here the 
initial tube pressure is higher at a given non-dimensional thrust. It 
is evident that high values of the non-dimensional thrust can be 
achieved before dissociation becomes important and that increasing the 
initial tube pressure P-| decreases the amount of dissociation in the 
propel lant. 

Figure 6 shows the effect of non-dimensional thrust on the fraction 
of the total thrust that is produced by the differential pressure-area 
product across the ramjet. The conventional ramjet operating with a 
perfect gas derives a smaller fraction of thrust from the pressure-area 
differential than the other ramjet configurations. The electrothermal 
plug ramjet operating with a perfect gas obtains a greater fraction of 
thrust from the differential of the pressure area product than the con- 
ventional ramjet for a given value of the non-dimensional thrust. The 
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electrothermal ramjet analyses which consider non-perfect gas effects 
show that including these effects causes the pressure induced thrust 
fraction to rise over that given for the perfect gas, electrothermal 
ramjet. The high tube pressure, high thrust ramjet operates with a 
lower pressure-induced thrust fraction than the low tube pressure, low 
thrust ramjet. 

Figure 7 illustrates the relationship between the exit/inlet 
static pressure ratio and the non-dimensional thrust for the different 
ramjet models considered. By design, the exit/inlet pressure ratio is 
unity for the conventional ramjet. The electrothermal plug ramjet 
operating with a perfect gas exhibits lower pressure ratios than the 
electrothermal plug ramjet operating with a reacting propellant. For 
the electrothermal plug ramjet operating with a reacting propellant, 
high thrust, high pressure operation is seen to exhibit slightly lower 
exit/inlet pressure ratios than low thrust, low pressure operation. 

The exit/inlet velocity ratios for the different ramjet configura- 
tions are shown in Fig. 8. The conventional ramjet exhibits a nearly 
linear dependence between the exit/inlet velocity ratio and the non- 
dimensional thrust. It has been suggested previously that, for a con- 
ventional ramjet, both the thrust (Eq. 9) and the exit/inlet velocity 
ratio (Eq. 8) show rates of change approximately proportional to the 

square root of the total temperature ratio (Tt /T- ), so a linear depend- 

'7 h 

ency between the thrust and the exit/inlet velocity ratio is to be ex- 
pected. The results obtained for the electrothermal plug ramjet 
operating with a perfect gas show that the velocity ratio is degraded 
with respect to the conventional ramjet. This curve is not linear 
because the outlet/inlet static pressure ratio increases with 



EXIT/INLET PRESSURE RATIO R,/P, 


V, = 7000 m/s 


20r T, = 300 °K 
S = 30® 

= 0.805 
d^/d, = 0.403 


ELECTROTHERMAL PLUG RAMJET ( A 7 /A,= I ) 


WITH DISSOCIATION 

= 9.5 X 10® Pa NO DISSOCIATION 




yy 


V% 


CONVENTIONAL RAMJET ( P 7 ./P, = I ) 
NO DISSOCIATION \ 


NON-DIMENSIONAL THRUST [ F / P, A, ] 

Fig. 7 Outlet/Inlet Pressure Ratio Profiles 


OF POOR 


: 1 : ' \X *-.v * ‘ * ' ' "' ‘ ■ ■'" ' V «' * 


i A 



Fig. 8 Outlet/Inlet Velocity Ratio Profiles 


3C 




KV 





t 


/ 


y“ 

f- 

!»< 


I. 


■■ ■ A- • dl 


non-dimensional thrust for the plug ramjet (a nearly constant static 
pressure ratio is essential to ensure linearity, as described by 
Eqs. 4 and 5). The electrothermal plug ramjet results including dis- 
sociation demonstrate that the exit/inlet velocity ratio is degraded 
with respect to the perfect gas configurations. Only one curve is 
shown for the cases including dissociation. This is because the data 
for the two different cases were very close, and a true distinction 
was difficult to show graphically. The actual data indicate that the 
high thrust high pressure device operates with a slightly greater velo- 
city ratio than the low thrust, low pressure device. This degradation 
in the momentum change of the propellant between the inlet and outlet 
is offset by greater exit/inlet pressure ratios so that desired thrust 
levels can be achieved. Figure 8 also brings to light the fact that 
the exit/inlet velocity ratio can be less than unity. This implies a 
momentum change of the propellant detracts from the thrust and requires 
that the differential pressure-area product must supply more thrust. 

In Fig. 6, these regions of operation where the pressure-induced thrust 
fraction is greater than unity are observed at the low normalized 
thrusts where the velocity ratio V^/V^ is less than unity. These regions 
also correspond to the low efficiency regions illustrated in Fig. 4. 

In Fig. 9 the relationship between the exit/inlet temperature ratio 
and the non-dimensional thrust is illustrated for the different ramjet 
configurations. The conventional ramjet exhibits a lower temperature 
ratio than the other ramjet configurations because of the more complete 
nozzle expansion. In the case of the electrothermal plug ramjet, the 
dissociating propellant cases are shown initially to exhibit greater 
temperature ratios than the electrothermal plug ramjet operating with a 
perfect gas, and the temperature ratios for the dissociating propellant 
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drop below ttioseshown for the perfect gas electrothermal plug ramjet 
as the non-dimensional thrust is increased. This lowering of the exit/ 
inlet temperature ratio is due to the dissociation of the propellant. 

The effect of changes in the flight velocity on the instantaneous 
efficiency can be seen by comparing the results of Fig. 4 (V-j = 7,000 m/s) 
with those of Figs. 10 (V^ = 15,000 m/sec) and Fig. 11 (V^ = 4,000 m/sec). 
All of the other operational parameters mentioned previously, such as 
the diffuser cone half-angle, are the same as for the 7,000 m/s flight 
velocity, except for the diameter ratios d^/d-j and d^/d^. The diam- 
eter ratio cl^/d^ = 2 has been specified to be constant for any velocity. 

As the ramjet flight velocity increases, the conical shock angle will 
decrease if the conical diffuser half-angle is held constant. Because 
it has previously been specified that the intersection of the oblique 
shock and the tube wall occurs at the diffuser shoulder, the tube diam- 
eter must decrease as the projectile velocity increases for a constant 
ramjet annular diffuser geometry. This is why the diamete*" ratios d,/d-| 
and d^/d^ increase as the flight velocity increases. Examining Figs. 10 
and 11 one sees that the same trends associated with the 7000 m/s 
results are evident in these figures (i.e., there is a value of the 
non-dimensional thrust for which the instantaneous efficiency is 
optimized). A very important feature of ramjet operation is brought 
out by examining Figs. 4, 10 and 11. This feature is that in order for 
ramjet operation to remain at a reasonable efficiency level during ac- 
celeration, the non-dimensional thrust must vary greatly over the 
region of acceleration. Assuming that a constant thrust level (F) is 
desired and that the inlet area (A^) does not change appreciably during 
the launch of a given plug projectile makes it apparent that the static 
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propellant pressure itmiediately upstream of the ramjet must decrease as 
the ramjet velocity increases so that efficient operation can be main- 
tained. It is noted that the data in Fig. 10 for the dissociating pro- 
pellant, electrothermal plug ramjet are terminated at a value of 90 for 
the non-dimensional thrust because ionization of the propellant started 
to become significant past this value. 

The results of Fig. 12 were obtained from a series of plots like 
those in Figs. 4, 10 and 11. They illustrate the variation in the non- 
dimensional thrust which must be achieved to attain the maximum instan- 
taneous efficiency over a velocity range of 4,000 m/s to 15,000 m/s. 

The non-dimensional thrust is observed to vary from a minimum of 4 at 

4.000 m/s to a maximum of around 50 for a velocity of 15,000 m/s. 

These results are for the dissociating propellant. 

Figure 13 illustrates the maximum instantaneous efficiency attain- 
able for a mission where a projectile is accelerated from 4,000 m/s to 

15.000 m/s utilizing two different thrust levels. The operational par- 
ameters are the same as previously defined. To attain the performance 
shown in Figure 13 the non-dimensional thrust must follow the non- 
dimensional ized thrust profile given in Fig. 12. This means that 

must be varied continuously along the launch tube. For high velocities, 
the high thrust profile exhibits performance superior to the low thrust 
profile because high static tube pressures are associated with high 
thrust operation, and propellant dissociation is less significant for 
the high pressure case than for the low pressure case. The optimum 
value for the instantaneous efficiency is observed to be near 7,000 m/s 
because this is the velocity for which this particular diffuser 
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performs the best. Changing the diffuser geometry will change the 
location of the optimum point. 

Several profiles are shown in Fig. 14 which detail the performance 
of a ramjet as it accelerates 10 kg at a level of 30,000 g's. These 
profiles were computed for various conical diffuser half-angles and for 
tube static pressures that remain at fixed values along the tube length. 
For purposes of comparison, the optimum efficiency performance data for 
the high thrust per unit area values given earlier are plotted 
with the results obtained for constant tube pressure acceleration. 

These values of the thrust per unit area coupled with the specified pro- 
jectile mass produces an acceleration of 30,000 g's. The early portions 
of the curves where no data are shown are those portions of the launches 
where the ramjet propulsion system could not achieve the specified 
thrust and where light gas gun induced acceleration would be used. The 
influence of the tube pressure can be clearly seen by comparing the re- 
sults for tube pressures (P^ of 30 and 60 atmospheres and for a conical 
half-angle of 30°. Early in the Operation, a high tube pressure is 
shown to exhibit superior performance, while towards the end of the ac- 
celeration period the lower tube pressure exhibits the superior per- 
formance. The results plotted for the optimum performance suggest that 
the variation in tube pressure is important in order to maintain high 
efficiency. 

The effect of varying the conical diffuser half-angle is also il- 
lustrated in Fig. 14. Having a diffuser cone half-angle that is either 
too shallow or too steep degrades the performance. For the missions 
analyzed in this paper, the overall optimum conical half-angle was 
near 30°. 
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By applying the first law of thermodynamics and the momentum 
equation to the control volume shown in Fig. 3b, expressions for the 
non-dimensional thrust and the instantaneous efficiency can be put in 
the following forms respectively 


P _ PiVjCV^/V, - 1) + Pi(P7/Pi - 1) 


FV^ p^V^CV^/V^ - 1) + P^(P7/P^ - 1) 
Q 


( 12 ) 


03 ) 


Examination of Eqs. 12 and 13 illustrates that for given constant inlet 
conditions the non-dimensional thrust and the instantaneous efficiency 
will remain constant if the heat addition per unit mass (q), the exit/ 
inlet velocity ratio (V^/V^) and the exit/inlet pressure ratio (P 7 /P-|) 
remain constant. In order for these parameters to remain constant it 
is required that the propellant thermodynamic properties remain con- 
stant throughout the ramjet flowfield shown in Fig. 3b. To ensure a 
consistent solution for the flowfield at state 4 shown in Fig. 3b, it 
is necessary to keep the diameter ratio d^/d^ and the conical diffuser 
half-angle (6) constant for whatever flight velocity is being consider- 
ed. Keeping these variables constant allows the results given here to 
be scaled for different size projectiles and thrust levels. 



CONCLUSIONS 


The results presented in this thesis have extended the results 
obtained from earlier analyses.^ The earlier analyses assumed isen- 
tropic diffuser operation and perfect gas propellant behavior, while 
the present analysis takes into account such effects as diffuser shock 
losses, dissociation and ionization of hydrogen propellant, and non- 
equilibrium processes in the supersonic nozzle expansion. It has been 
demonstrated that diffuser losses, while significant, do not degrade 
ramjet performance to unacceptably low levels. The diffuser geometry 
does affect performance significantly if the diffuser conical half- 
angle is too shallow or too steep. At the overall optimum conical 
half-angle for the cases presented here (''-30°) however the performance 
degradation is not excessive. 

In the analyses conducted in this thesis, it has been assumed 
that the propellant is in thermodynamic equilibrium up to the throat of 
the annular nozzle, and nonequilibrium processes are allowed to occur 
in the supersonic nozzle expansion downstream of the nozzle throat. 
Chemical nonequilibrium in the nozzle expansion was evaluated, and it 
was found that the nozzle flowfield was typically found to be very near 
chemical equilibrium. Because chemical relaxation times are generally 
larger than other nonequilibrium process relaxation times (such as 
vibrational, translational and rotational processes) the assumption 
that the propellant is in thermodynamic equilibrium should not intro- 
duce significant errors into the analyses. 
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It has been demonstrated that the static propellant pressure im- 
mediately upstream of the electrothermal ramjet affects performance 
significantly. Changing this pressure effectively changes the mass 
flow rate and the rate of heat addition per unit mass and it is demon- 
strated that for a particular flight velocity there is an optimum rate 
of heat addition per unit mass. When dissociation and ionization are 
allowed to occur, it is shown that decreasing the static propellant 
pressure increases the losses due to dissociation and ionization. 

However, the results presented demonstrate that it is feasible to oper- 
ate in operational regimes where dissociation and ionization do not 
significantly degrade ramjet performance. To maintain peak efficiency 
the static propellant pressure immediately ahead of the ramjet must de- 
crease as ramjet velocity increases. 

The assumption of quasi-steady fluid flow used in this analysis 
will surely introduce some error. It is possible that the mass of the 
propellant contained within the ramjet could change, and the fact that 
the propellant within the ramjet is also accelerating could also in- 
troduce significant errors. The acceleration could cause a "hydrostatic" 
pressure difference to occur within the ramjet heat addition zone, 
which would possibly invalidate the results presented. An initial 
analysis of these effects has been conducted and is detailed in Appen- 
dix r. This analysis demonstrates that the non-steady gas dynamic 
effects just mentioned do not introduce errors which invalidate the 
results which are presented. The results from this analysis are based 
on calculations which compute flowfield changes over typical ramjet 
fluid particle residence times. Indeed, the idea is put forv/ard that 
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proper tailoring of the upstream propellant conditions will minimize 
the unsteady gas dynamics effects to the point where they are not 
important. 

Much analysis remains to be done to investigate completely the 
effects of all of the parameters of concern. The problems of projectile 
material selection and design have not been addressed nor have the de- 
tails of the heat addition process been examined. It is possible that 
pressure oscillations could occur in a device of this type, and these 
oscillations will be related to the unsteady gas dynamics effects. 
Additional work related to maintenance of the correct projectile base 
pressure and static pressure immediately ahead of the projectile remains 
to be done. Although this aspect of the problem has not yet been 
studied, it intuitively seems possible that methods will be found which 
will allow for maintenance of the correct pressures. It is also possible 
that three-dimensional flow effects could become significant. 

The results presented earlier have demonstrated that the conven- 
tional ramjet has the capability to outperform the plug ramjet. However, 
the mechanical complexities associated with conventional ramjet operation 
(such as providing efficient shock swallowing and producing the necessary 
nozzle conditions) are very significant. It is therefore argued that 
the plug ramjet, because of its apparent mechanical simplicity, is a 
more desirable device for effecting high acceleration launches. 

The results presented suggest that the electrothermal ramjet con- 
cept cannot be rejected on the basis that it violates a basic principle 
of fluid mechanics or thermodynamics. Further, the results presented 
suggest that the electrothermal ramjet has the potential of acceler- 
ating payloads to high velocity with a good energy utilization 
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efficiency. For these reasons It is concluded that further analysis 
related to the design and development of the electrothermal ramjet 
concept is warranted. 
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APPENDIX A NOMENCLATURE 

A = area, m^- 

Cp = constant pressure specific heat, J/Kg°K 

h = enthalpy, J/Kg 

m = mass flow rate. Kg/s 

n = normal to a control surface 

P = static pressure. Pa 

P^ = stagnation pressure. Pa 

R = radial distance, m 

r = radial direction subscript 

T = static temperature, ®K 

Ty = stagnation temperature, ®K 

V = velocity, m/s 

6 = conical diffuser half-angle, radians 

Y = ratio of specific heats 
p = density, Kg/m^ 

Oj = conical shock half-angle, radians 

A bar over any variable indicates that this variable represents an 
average value, a numerical subscript denotes a position shown in 
Figure A1 , and a subscript of r or e indicates a direction as specified 
in Figure A2. 





APPENDIX A 


CONICAL FLOWFIELD AVERAGING 

In this analysis it is assumed that the working fluid behaves as a 
perfect gas, and that the flowfield is axis-symmetric and quasi-steady. 
Assuming that the working fluid behaves as a perfect gas should not 
introduce too much error because for shallow cone angles the flowfield 
conditions should not change significantly across a conical shock even 
for high incoming Mach numbers. The quasi-steady assumption is justi- 
fied in Appendix D where it is shown that fluid particle residence times 
are small compared to the times required for flowfield changes to be 
significant. The geometric constraints should not let the ramjet pro- 
jectile deviate significantly from axis-symmetric conditions, so 
assuming axis-symmetric flowfields exist should not introduce signifi- 
cant errors. 

Consider a conical diffuser of the type shown in Fig. A1 to be 

operating at a zero angle of attack within a launch tube. Assume that 

the projectile is traveling at a high enough velocity so that conical 

shock attachment is assured. The conical shock formed will intersect 

the tube wall as shown in Figure A1 . For a given conical half-angle (6) 

the flow from 1 to 2 can be calculated given conditions at 1 using the 

A1 

Taylor-Maccol method. Indeed, the conditions anywhere within the 
axis-synmetric flowfield can be calculated using this method. The 
problem is to describe the flow between 1 and 4 for given conditions 
at 1 and a specified diffuser geometry. In reality a pattern of re- 
flected oblique shock waves would exist downstream of the initial 
conical shock depending on the annular diffuser throat area and the 
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pressure at 4. If it happens that the flow is restricted at the dif- 
fuser throat or downstream of 4 then a bow shock could exist ahead of 
the diffuser. If the downstream conditions are selected as they are 
for designs considered here so that the flow is not restricted and that 
the back pressure conditions are correct then a reflected shock patte'^n 
terminated by a normal shock will exist. The location of the termina- 
ting normal shock can be dictated by varying the diffuser throat area 
and the pressure at 4. 

For hypersonic flows, because of separation and boundary layer- 

shock wave interaction effects, diffusers perform close to what is 

A2 

predicted for normal shock diffusers. Thus in this analysis the 
following assumptions are used: 

1) T.'ie separation and boundary layer-shock interaction behavior 
does not become important until the reflected shock pattern 
begins downstream of the intersection of the tube wall with 
the initial conical shock. Neglecting separation and boundary 
layer effects in this initial conical flowfield region should 
not produce large errors because the scale of the diffuser is 
expected to be small. Neglecting these effects also allows 
the Taylor-Maccol method of analyzing conical flowfields to 

be utilized. 

2) The flow downstream of the initial conical flowfield is 
rapidly mixed and can be modeled by a one-dimensional flow 
which satisfies the conservation laws. It is noted that the 
conservation laws are insufficient to completely specify the 
one-dimensional flow; an additional assumption which will be 
described shortly is also needed. 



3) A normal shock terminates the one-dimensional supersonic 
flowfield. 

When these assumptions are applied the flow processes in the dif- 
fuser can be modelled in the following way suggested by Fig. A1 . The 
free stream flow passes from state 1 through the conical shock into the 
conical flowfield of region 2. The properties of this flowfield, which 
vary from point to point and must be calculated using a numerical pro- 
cedure, are averaged over the surface immediately upstream of the normal 
shock to obtain the properties at 3. These averaged properties serve 
as input to the one-dimensional normal shock relations and yield the 
properties of the fluid immediately downstream of the normal shock at 
state 4. 

The additional statement needed to connect the conical and one- 
dimensional flowfields will now be considered. The choice used here 
is made through the definition of a mass flow rate averaged velocity 
for the conical flowfield and equating this to a characteristic one- 
dimensional velocity V3. In Fig. A1 let states 2 and 3 represent the 
conical and one-dimensional flowfield states respectively. The mass 
flow rate average velocity at 3 is defined by 

r J |Vl P(V • n) dA 

V3 = 

c^s/pCV • n) dA 

where the subscript c.s. refers to any arbitrary control surface which 
the complete conical flowfield passes through. The integration is 
simplified if a spherical coordinate system is placed at the vertex of 
the physical cone and the control surface is taken to be a spherical cap 



as Is shown in Figure A2. 
For this choice 
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V * n = = velocity in the radial direction 

|V| ■ lv‘ * 

p = p(o) [for conical flow] 
dA = 27 tR 2 Sino de . 

Thus it follows that 

r P (V^ + Vy 2)’-5 27tR Sine de 

p 27 tR 2 Sine do 

6 

In this expression the factor 2 tiR 2 appears in both the numerator and 
denominator and thus will cancel. This implies that the one-dimensional 
properties of a conical flowfield are independent of the scale of the 
flow conditions. It follows then that the value for the average velocity 
immediately upstream of the normal shock is calculated from 

js p V^(V2 + V2)’^ Sine de 

. 

jS p Sino do 

6 

Hass continuity, which must also be satisfied in the conical flowfield, 
requires that 

_ _ Os 

P 3 V 3 A 3 = / 2 itR 2 p Sino do = m 3 
6 
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where 

6 - 

A3 = 2 tiR 2 / Sine de = - 2 itR 2 (Coseg - Cos6). 

6 

Combining these expressions one obtains an expression for the mean 

density immediately upstream of the normal shock. 

e. 


M3 


/ p V Sine do 
6 


V3 (Cos6 - C0S65) 

The first law of thermodynamics must also be satisfied. 


^T3 '^Ti “ S ^3 ^ ^ 3 ] 


Solving for the mean temperature upstream of the nomal shock yields 

hx 


rr2 

xrJ, „ 
y "2R 


The ideal gas equation is assumed to be valid. 




P3 - P3 R T3 

The average stagnation pressure is defined by 



These average values for the flowfield variables are used as input for 
the normal shock relations. 

The only criterion applied to determine the annular diffuser throat 
area is that this area be large enough to swallow the flow downstream 
of the normal shock. 



The Taylor-Maccol method of analyzing conical flowfields, which 
was utilized in the analyses conducted in this thesis* is not presented. 
For a complete description of this procedure see reference Al. 
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APPENDIX B 

EQUILIBRIUM HYDROGEN PROPERTIES 





APPENDIX B NOMENCLATURE 


A = Helmholtz function, J 

= diatomic hydrogen dissociation energy, 7.171 x 10"^® 

= lowering of the ionization potential for specie i, J 
i . 

e = electron charge, 1.6021 x 10"^^ Coulomb 
g = degeneracy 

h = Planck's constant, 6.62517 x lO"^** J-sec 
h = enthalpy per unit mass, J/Kg 

h = enthalpy per unit mole, J/Kmole 

K = Boltzmann's constant, 1.3804 x lO"^^ J/°K 
M^. = molecular weight of specie i, Kg/Kmole 

m^ = mass of specie i 

N. = number of particles of specie i 
n^Q^ = number density of ions, m"^ 

P = pressure, Pa 

Q^- = charge of specie i 

R = universal gas constant, 8.314 J/mole°K 

S = entropy, J/°K 

s = entropy per unit mass, J/Kg°K 

s = entropy per unit mole, J/mole°K 

T = temperature, °K 

U = internal energy, J 

u = internal energy per unit mass, J/Kg 

u = internal energy per unit mole, J/mole 

V = volume, m^ 

v^. = molar specific volume of specie i, m^/mole 




= mole fraction of specie i 
= partition function for specie i 

X 1- 

= electronic energy of the i level for a given specie, J 
= electron affinity for the i^^ specie, J 
e^. = ionization energy for the i specie, J 




APPENDIX B 


EQUILIBRIUM HYDROGEN PROPERTIES 


In analyzing ranyet operation it is necessary to calculate the 
equilibrium composition and thermodynamic properties for hydrogen. This 
is done by the methods presented in this Appendix. The species con- 
sidered in this analysis are 


Ha. H. H~. e" . 


When a mixture of these species is at equilibrium the Helmholtz 


function 


A = U - TS 


is at a minimum. In this analysis it is assumed that the Helmholtz 
function can be expressed in the form 


* ■ "pg * "c 


Where Ap^ is the perfect gas Helmholtz function, which is valid for 
noninteracting particles, and is the coulombic interaction Helmholtz 
function which comes into play for the case where the mixture is ionized. 
The coulombic term can be considered as a correction to the perfect gas 
form. The perfect gas term can be expressed in the fonn 


Apg = - KT [E N.(in(-^-^) + 1)] 
’ i 


where N. is the number of particles of the specie and (?^/N.)is the 
perfect gas partition function for the i^^ specie. The coulombic term 


i 
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can be expressed in the form 




where is the charge of the i'^” specie and e is the charge on an 
electron. At equilibrium 

dA = dApg + dA^ = 0 . 

The differentiation is taken to occur at a constant temperature and 
volume. Variations in the Helmholtz function due to changes in the 
composition are considered. The forms for the differentials are 


dApg = - KT Z tn (Z^/N^) dN. 


and 


dA^ = -e3 (-^-2^) [iiN.Q^] [IQ? dN.] 


Combining these expressions yields 

This expression can be put into a slightly different form 


I tn dN^ + I 


I, 


-^(f) D:N.Q?] Q|dN. 


= 0 . 


Factoring out dN^ yields 


Z, 


£n 


1 . e (j) Q, 


i (KT) 


dN^ = 0 


or 


I J.dN. = 0 
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The chemical reactions considered are 

Pg H 2 ^ 2H a) 

Pb Hz - Hz"' + e" b) 

K ^ + e" c) 

p^ H^H + e . d) 

where the degree of reaction for each chemical reaction is defined as 
being the variable p. If each reaction is assumed to proceed slightly 
from left to right then the change in each component is given by 


^^Hz " ■^''a ■ % 
dNHz+ " ''"b 

d^H = 2 dpg -dp^ + dpj 
dN/ = dp^ 
dNn" = -dp^ 

dNe = dpj^ + dp^ + dp^ . 

By combining these expression, with the expression 
Z dN^ = 0 

the following equation can be obtained 

^ 2J^) dp^ + (-J^^ + Og.) dp(^ + 

+ dg-) dp^ + + dg-) dpj = 0 . 


Since the cnemical reactions and hence the dp values are assumed to be 
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independent, each quantity in brackets must be equal to zero at 
equilibrium. This equation is used to help determine the chemical 
equilibrium composition. 

Consider the rccction 

H 2 2H . 

For this reaction the relation 
is used. In more complete form this equation becomes 



Zm h 2 

= 2.n/ + Arr-(J) [^NQ.] Q 
Nh (kT) ^ 1 1 M 

but 

%2 


so this equation degenerates to 


\ "«2 

which is the correct form for a conventional perfect gas analysis. 
Multiplying each side of the equation by 

N 

2 — 2 
Un.) N 


allows this equation to be put into the form 



v : 
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*1 

''«2 (Z/N)„^ 

where N is total number of particles, 
For the reaction 

H2 ^ H2^ + e 
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the relation 


-Jh 2 " 'H2" " ^e- = 0 


is valid. In its more complete form this equation becomes 

Z., 


H?. 

,n + - 


H2 (KT) 


777 (J) [ 2 :N^Q?] Qu = 


Zu + 


fn 


^e- 


^ trr ^ -”^7^ ( J ) [^n.q2] 


(KT) 




'I 




where 


%2 = ° 

^H 2 " = ^ 

Qe = - 1 . 

After rearrangement this equation degenerates to 



= 2 e 


.^5 


(KT) 


n ^ 

T 7 ? "ion 


"ion 


IN.Q? 


fr" 


V- 
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and represents the ion number density. Rearranging and multiplying 

through by N/N^ yields 


^Hre- 


(Z/N)„+ (Z/N)g- 


exp 


(KT)'* 



For the case where n^^^ is a small number, the exponential term on the 
right hand side of the above equation will be close to one, which cor- 
responds to the form obtained when only a perfect gas analysis is done 
(i.e. the Saha Equation). When the number density of ions is signifi- 
cant, however, a coulombic field will be present, and this field will 
affect the ionization of the different particles. The net effect of 
this coulombic field is to decrease the amount of energy needed to re- 
move an electron from the influence of a given nucleus (the ionization 
energy is decreased). Consequently, the exponential term in the 
preceding equation can be thought of as representing an expression for 
the lowering of the ionization potential. The lowering of the ioniza- 
tion potential can therefore be expressed as 


aE. = 2e 
ion 

H, 


3 *^ion 


KT 


and the expression for the mole fraction ratio becomes 

'AE 


Y + Yg- 

H2 


(Z/N)H^(Z/N)g- 


(Z/N) 


expl 




ion 

H2 

KT 


For the chemical reaction 
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the same form is obtained for the expression for the mole fraction 
ratio namely 


Yu+ Y - (Z/N)u+ (Z/N)^. 


H__e 


H'' 
(Z/N), 


exp 


''A E. 

ion 


KT 


where 


^^ion - 2 eM^) 


H 


H, 


An interesting result is obtained when the reaction 


H' 4 H + e" 


is considered. For this reaction the relation 




is valid. When this equation is expanded it assumes the following 
form 


where 




2nj^ +J!,n|^ + ~TT2 ii^) [^^N-jQ-j] Qe~ 

^e- (KT) ' ^ lie 


Q^_ =Qfi-=l . 


The coulombic interaction terms will cancel so this equation assumes 


Z, Z^- _ N, W,. 




N,- 


the form 
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which predicts that, for a negative ion, there is no lowering of the 
ionization potential. However, if the negative ion was doubly charged 
and only one electron was removed, then a change in the ionization po- 
tential would be predicted. Multiplying the previous equation by N/N^ 
yields 

Y^Y^. (Z/N)^^ (Z/N)^- 


'H' 


(Z/N). 


At this point in the analysis, four equations describing the rela- 
tionships of the mole fraction species have been derived. In summary, 
these equations are 




H 


(Z/N). 


Ynt Ye- (Z/N)h. (Z/N)^ 


(Z/N), 


exp 


^ ^ion 
Hz 

KT 


V V _ 


(Z/N). 


exp 


A E 


ion 

H 


KT 


Yh Y^- _ (Z/N),^ (Z/N), 




(Z/N)^- 


However, because there are six species considered, two more equations 
must be obtained in order to ensure a unique solution. One of these 
two equations can be derived by assuming that charge is conserved 


’'h" * V2 - ’'h- = » • 
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The second equation comes from summing up the mole fractions 


V*V V'V-^V 



The variable (Z/N)^ is the partition function for the i^" specie 
evaluated at the total system pressure and temperature. In order to 
evaluate the partition function it is necessary to select an energy 
datum, and the datum selected for this analysis is the ground electronic 
state of the H 2 molecule at 0°K. The- partition function of each specie 
is computed relative to this state. 

The partition function for any specie i can be expressed in the 



where the translational contribution is given by 

3/2 

(KT) 


(V"), 



V2 


The form of the internal contribution varies depending on whether the 
specie under consideration is diatomic or monatomic. A diatomic specie 
will have rotational and vibrational contributions, and the diatomic 
species and some of the monatomic species will have electronic contri- 
butions. The possibility of having a chemical energy contribution to 
the internal partition function also exists for some of the monatomic 
and diatomic species, notably H 2 > H, H , H', and e", depending on 
where the ground state energy of the given specie relative to the 
datum selected above lies. 



r 





Consider the monatomic specie of neutral hydrogen (H). The in- 
ternal partition function for H is given by 



7 =7 7 

internal chem 

elec 


H H 

H 


/ Dq 

\ 

where 

^chem " (■ 2KT 

) 


H 

/ '^e 

and 

^elec ~ ^ ^e. 
H J 

f--i 

\ KT 



This expression assumes that each H atom receives 1/2 of the chemical 
energy required to dissociate a ground state H 2 molecule. The de- 
generacy (g„ ) in the electronic contribution relation is given by 

2 

9g = 2j j = 1, 2, 3, . . . 


and the energy (e ) is given by 

j 




h2 


j = 1 , 2 , 3 , . . . 


I 


In this equation for the energy state, L is the charge of the nucleus 
(1 for H), e is the electron charge, m is the electron mass, h is 
Planck's constant and C is a proportionality constant. Inserting these 
relations into the expression for the electronic partition function for 
the hydrogen (H) atom yields 


^elec 




27(^C^L^e'*m -i 42 


h^KT 



f 


Examination of this equation brings to light the fact that the series 

will diverge unless a cutoff (J) for the summation is specified. A 

method of calculating this cutoff will be given later in this Appendix. 

+ • + 

The other monatomic species are H , e and H . For the H particle 
the value for the internal partition function is 1, and the chemical 
partition function is given by 


Wrn “ '’'P 


- (D /2 + c, ) 
0 1 ^ 


In this expression e- is the ionization energy of a hydrogen atom. 

+ » 

The H particle is assumed to carry all of the ionization energy. For 
the electron (e") the internal partition function is 


because of the spin contribution to the degeneracy. The electron is 
assumed to carry no chemical energy so the chemical contribution to the 


partition function is 


= 1 . 


The H' ion has an electron affinity value that must be taken 

into account when calculating the chemical partition function. The 
chemical partition function is given by 


H" 


" Vfin) . 


Only the ground electronic state is considered for the H“ ion and 
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the internal partition function (electronic) is 

Zi = 1 . 

int 

H" 

Higher electronic states are neglected because the existence of elec- 

B 1 

tronically excited negative ions is rare. The binding energy for a 
negative ion is generally small and it is not too difficult to renxive 
an electron from a negatively charged particle. Also, because the 
binding energy is small, it seems to make sense to neglect electronic 
states above the ground state because the higher electronic states ef- 
fectively reduce the amount of energy required to remove an electron. 

The binding energy for the extra electron on a negative ion is not large 
because the inner electron(s) repel the extra electron with a charge 
that is equal to the charge of the nucleus, effectively shielding the 
extra electron from the nucleus. 

For the diatomic species (H2, H2"^) there are rotational and vibra- 
tional contributions to the internal partition functions that must be 
taken into account. In this analysis coupling between the rotational 
and vibrational en'‘rgy storage modes is allowed, and an expression 
characterizing the rotational and vibrational modes for separate elec- 
tronic levels is presented later in this Appendix. The complete in- 
ternal partition function for the H2 molecule (and H2'*') is expressed by 


^int -J, ^r» ^e. 

where the summation must be terminated at a finite electronic energy 
level. The internal partition function is more completely expressed in 



__ ,i^ j|^»Wl Wi W il| l |l 1. ll ^^ ■ 
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the form 


Z, “I 9g exp ( — A 

int an ®j \ KT / 


-1 (J_„) 


\ .- 


1 + y. + ^ + -A — , + t 

^ ^ -1 (e^-l )2 

In this expression o = 2 because a diatomic molecule is syrmietric. The 
following parameters are evaluated from the rotational and vibrational 
constants for the given electronic levels (the rotational and vibra- 
tional constants can be found in compilations such as the one in 
Reference B 2 ). 


he 


Cw - X w j 


KT -e e e' 


X W 

e e 


6 = 


Y = 


w 

- X W 

e 

e e 

he 


w 

(B^ - 1/2 


a 

B~ 

- 1/2 a 

6 

£ 


D 


£ 

B~ 

- 1/2 a 

e 

e 


X W , B 
e £ 1 


e e 

coiistants and is computed using 

4 B 


D = 






The H 2 '*' molecule has a chemical (ionization) energy that must be 
taken into account, and the chemical partition function for is 
given by 


^chem 

H2+ 


exp 


KT 



It should be explained here that the ionization energies used in eval- 
uating the partition functions are values published in a reference such 
as Reference B3. It is assumed that these values were obtained for 
rarified, very weakly ionized gases where coulombic interactions are 
unimportant. The calculation of the ionization energy for a given 
specie takes the published value as a baseline, and this value is then 
changed in the manner detailed earlier where the changes in ionization 
energies are defined for the different species. 

In the evaluation of the various partition functions it is neces- 
sary to specify cutoff values whenever electronic levels other than the 
ground level are considered. In essence, the problem is establishing a 
criterion for when an election can no longer be considered to be bound 
to the nucleus of an atom or to the nuclei of a molecule. The criterion 
selected for this analysis is to not allow any electronic energies which 
exceed the ionization energy (taking into account the lowering of the 
ionization energy discussed earlier). 

Solving for the equilibrium composition is an iterative process 
because the equations are nonlinear and the lowering of the ionization 
potential must be determined implicitly. Tiie method used in this 
analysis is to assume initially that no ions are present and to apply 
the model just presented and compute the ion density. A new ion density 
based on the computed value is then used as input and this process 
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continues until the difference between the assumed and computed ion 
density is sufficiently small. The value required for the number density 
of the ions is computed from 



_^on 


N ^ion 


ion 


_P 

Kf 


(Y^. * V„. * * Y„ ♦) 


Once the equilibrium composition has been found, the thermodynamic 
properties of the mixture must also be evaluated. It is noted that the 
coulombic corrections previously referred to in this Appendix will 
affect the composition and the thermodynamic properties of the species 
because the ionization energies and the allowed electronic excitation 
levels are changed by allowing coulombic corrections to appear. Also, 
the coulomb field itself will possess definite thermodynamic properties 
and these properties will contribute to the thermodynamic state of the 
complete hydrogen mixture. For example, the ions will be slightly 
attracted to each other so that the actual system pressure will be 
slightly less than is predicted by the ideal gas equation. However, in 
the spirit of Reference B4, it is assumed in this analysis that the con- 
tribution of the coulombic field to the total system themwdynamic state 
of a hydrogen mixture is small and can therefore be neglected without 
introducing significant error. 

The entropy (?), enthalpy (h) and internal energy (u) of a mixture 
can be expressed on a molar basis by the following relations 


s" • J V, 5, 
h - I Y, h, 
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Dividing each equation through by the average molecular weight which is 
defined by 

“ y V. M. 

allows these properties to be expressed on a per unit mass basis 

Vj 

' ■ i '’'r«7 
h - II' A 

): "i "i 

Ih 

The internal energy on a molar basis of a specie i (u^ can be 


expressed in the form 


3 tn Z. 


u. = RT2| - 


3 T 


Since the enthalpy (h^) is defined by 


h, . Uj + p. V, 


where v^ is the molar specific volume, and the ideal gas equation 




is assumed to be valid, the enthalpy per mole can be expressed in the 
form 

, / 3 an Zi 


+ RT . 
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The entropy per unit mole is expressed by POOR QUALITY 


s, = R 


h I “i 

n7 ^ 1 I ^ r 


The variable R is the universal gas constant. 

In order that the thermodynamic properties can be calculated the 
partial derivative appearing in the thermodynamic relations is expressed 
in the form 

D (an 

TT 


(IL-V) 


1 f ^ 


Z. V 8 T 

V-. ’ ' 'V-. 

It is possible to separate the translational partition function out of 
this expression. The partial derivative of the translational partition 
function for any specie is given by 


1 I ^ ^trans 


*trans' 


3 T 


- 3 1 

- 2 j 


^i 


Consequently, the translational contributions to the internal energy and 
enthalpy are given by 


u. = 3/2 RT 
trans 


h. = 5/2 RT 


trans 

When the chemical energy is considered, the internal energy and enthalpy 
are given by 

“i =^i 
chem chem 


wh 


lere is the dissociation energy and is the ionization energy. 
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The electronic contribution to the themiodynami c preoperties must 
be accounted for in the monatomic hydrogen (H) atom. The form for the 
partial r'erivative is 


nj2-l)exp 

u — = -L- J =J ^ 

3 T KT^ J / 

I exp( 

j=l ' 



The value obtained in this summation is utilized to calculate the 
enthalpy (h) and the internal energy (iT) for the hydrogen atom. The 
constant F is given by 


p = Zn^C^L^e^^m 
h2 


For the diatomic species the internal contributions to the parti- 
tion function are given by 

J 

Consequently the partial derivative is evaluated in the form 



For a given electronic energy e- the partial derivative of the elec- 

J 

tronic contribution function is given by 
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The expression for the partial derivative of the rotational -vibrational 
term is given below 


^ (^rv) 

~TT~ 


®i 

= AB + 


CD 


where the parameters A, B, C, and D are evaluated using 


A 


1_ 

ay 






4x*t^e^ 

(eTT)T 


C 




1 


X 2x*t 

(e^-1)' 


D 


1 l +e"^ (t- 1) 


Figures Bl, B2 and B3 illustrate results obtained for the mole fractions 
of the major constituents of equilibrium hydrogen for different pres- 
sures and temperatures. The constituents H 2 ’*’ and H” were never signifi- 
cant so results for these species are not shown. It is noted that even 
for the 10 atmosphere case (Fiq. Bl ) ionization does not begin to become 
significant until the temperature approaches 12,000 Kelvin. As a test, the 
calculated values for the thermodynamic properties and the composition were 



MOLE FRACTION 





compared to values obtained from Reference B5. In the cases evaluated, 
the percent differences between the results from Reference B5 and the 
results calculated were typical 1 of the order of 1% to 3%. It is recog- 
nized that the results of Reference B5 doe not include ionization. How- 
ever, the excellent agreement between the low-temperature results 
presented in Reference B5 and those given by the model described here 
is assumed to indicate high accuracy in the high- temperature regimes as 
well . 

The electrothermal ramjet was typically found to be operating with- 
in internal propellant tSmperature regimes of the order of 10^ Kelvin 
and pressure regimes of the order of 10^ atmospheres and it was found 
that ionization was typically unimportant for the ramjet operational 
conditions presented in this thesis. 
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APPENDIX C NOMENCLATURE 

B = constant used in Eq. 8, 78.01 
= mass fraction of H 

Do = H 2 dissociation energy, 7.171 x 10"i® J 

H = monatomic hydrogen specie 

H 2 = diatomic hydrogen specie 

K = Boltzmann's constant, 1.3804 x lO"^^ j/°k 

6 

Kj, = backward rate constant, 

m3 

= forward rate constant, 

M|^ = molecular weight of H, 

N = constant used in Equation 9 
Rj = dissociation rate, 

R = recombination rate, 

T = temperature, °K 

t = time, sec 

u = velocity, m/s 

V- = volume, m^ 

X = distance, m 

p = density, Kg/m^ 


harmonic oscillator constant for H 2 , 6338. 3°K 



APPENDIX C 


NONEQUILIBRIUM NOZZLE FLOW 


In the analysis of supersonic nozzles it is desirable to evaluate 

the nonequilbrium effects caused by finite rate reactions. This 

Cl 

Appendix describes the Bray sudden freezing technique which was 
utilized in the analyses presented in this thesis. 

The dominant species in this analysis are assumed to be H 2 and H 
and the following chemical reactions are considered 


H 2 + H 2 2H + H 2 (1) 

H 2 + H t 2H + H . (2) 

Four body and higher collisions are neglected. The rate equations 
which are assumed to apply to the above chemical reactions are for Eq. 1 


and for Eq. 2 


= 2Kf [H2]2 - 2K^ [H]2[H2] 
H2 H2 


(3) 


"Dt^ = 2K^ [H] [H 2 ] - 2K^ [H]3 . (4) 

H H 


[H] and [H 2 ] represent the molar densities of H and H 2 and K^andKj^ are 

Hz Hz 

the corresponding forward and backward rate constants pertaining to 
Eq. 1. Similarly K, and K. correspond to Eq. 2. Combining these 
equations to obtain the net molar production rate for atomic hydrogen 


yields 


Ml. ! 

Dt 1 

2Kf[Hz]2+ 2K^[H] [Hz] 

1 1 

i - i 
1 j 

2Kb[H]2[Hz] + 2Kj^ [H]3 

1 

Hz H 1 

1 1 

Hz H 


( 5 ) 
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Defining 

R. = 2K. [H^2+ 2K. [H] [Ha] 

H, » 

= 2K^^ [H]2[H2] + 2Kjj [H]3 
Ha ^ 

where represents the rate of dissociation and R^ represents the rate 
of recombination permits Eq. 5 to be expressed in the form 




fiOH = R _ R 
Dt "^d '^r 


( 6 ) 


It should be made clear at this point in the analysis that the previous 
equation assumes the volume of the hydrogen system does not change with 
time. The expression for the volumetric substantial derivative is 


.M ^ 

^ dt * 


The magnitude of this term is assumed to be small compared with R^ and 
and will therefore be neglected in this analysis. 

If the mass fraction of atomic hydrogen (C^^) is defined as the mass 
of H divided by the total system mass then the time rate of change of 
for a system can be defined by 


DC^ _ DLHJ 


= (R. 


R 1 ^ 

r' p 


(7) 


Dt Dt p ^"d 
where p is the system density (which is assumed to be constant) ind 
is the molecular weight of H. 

For a supersonic nozzle operating under quasi-steady conditions the 


substantial derivative expressing the time rate of change of the mass 
fraction of H is given by 
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DC^ _ 

"dF ' “ ~W 

where u is the local velocity and X is the distance downstream from the 
nozzle entrance. 

Consider a supersonic nozzle operating under quasi-steady condi- 
tions and assume that equilibrium flow exists up to the nozzle throat. 
Near the throat the local fluid velocities will be small, the chemical 
reaction rates will be large, and the characteristic chemical relaxation 
times will be small compared to the time that it takes a fluid element 
to move a fractional distance down the nozzle. Consequently the pro- 
pellant flowfield will be very close to chemical equilibrium, and for 
this analysis it will be assumed to be in complete equilibrium. As a 
fluid element accelerates down the nozzle, the molar densities of the 
species and the local static temperature, will drop and these two effects 
will affect tb..' flowfield significantly. The effect of dropping the 
molar densities cf the reacting species can be seen by examining Eq. 5, 
which shows that dropping the molar densities decreases the net reaction 
rate and this in turn increases the characteristic fluid particle 
chemical relaxation time. The effect of dropping the local static tem- 
perature is brought out by considering the forward and backward rate 
constants used in Eqs. 3 and 4. The backward rate constants used in 
Eqs. 3 and 4 are respectively 


_ 7.5 X 10’2 


Kmole"^ sec 


K. = 5.9 X 1Q13 

b„ — 


Kmole^ sec 
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For either Eq. 3 or Eq. 4 the forward and backward rate constants are 
related by the following expression 


-I = B 


exp - 


KT 


exp 


( 8 ) 


where 


B = 78.01 


Kii ole 


0 . 


51960 -X 


= 6338.3 '^K . 

This equation, which is actually a form of the equilibrium constant for 
the -► 2H reaction, can be derived by using the results of statistical 
mechanics and by assuming that electronic excitation is not important 
and that the H.-- molecule behaves as a harmonic oscillator. Examination 
of Eq. 8 illustrates that as the local static temperature drops, the 
ratio becomes small very quickly (this is due to the exponential 

term involving the chemical energy which dominates the expression). 
Ihis effectively causes to become smaller than (hence recombina- 
tion is favored). 

Using these results the Bray sudden freezing criterion can be 
stated in the following ways 

1) Near the throat the flow velocity is not large and the 
chemical relaxation rates are large because the local 
temperatures and the reacting specie molar densities are 
high. In this region the characteristic chemical relaxation 
times are small compared with the time that it takes a fluid 
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element to nwve a fractional distance down the nozzle and 
consequently the flow is assumed to be in chemical equili- 
brium. 

2) In the flow downstream of the throat, the molar densities 
and the temperatures drop while the flow velocity increases. 
The chemical molar densities lag farther and farther behind 
the equilibrium values and finite rate reactions begin to 
become important. 

3) The propellant eventually reaches a point where the net rate 
of chemical reaction is sufficiently small and the local 
velocity is sufficiently large so that the partial derivative 
DC^/3X becomes small. Downstream of this point in the flow- 
field the propellant is essentially chemically nonreacting 
and it is assumed to be fixed in composition (frozen). The 
characteristic chemical relaxation times are large compared 
to the time required to tranverse the flowfield. 

4) The transition from equilibrium to freezing occurs very 
quickly within a narrow region, and the location of this 
region is defined when the following relation is met 

.) C|i M,, 

„ -/ ^ R, » (9) 

where the subscript e denotes equilibrium conditions and N 
is a constant of the order of unity. This expression states 
that the narrow freezing region is located at the point where 
the local substantial derivative of the mass fraction of H 
is of the same order as the dissociation term appearing in 
Eq. 6. This implies that is beginning to become larger 



than Rj. Consequently, the flow is beginning to depart 
from chemical equilibrium and the characteristic chemical 
relaxation times are increasing. A value of unity was 
used for N in this analysis, which was suggested by 
Reference Cl . 

In the analyses conducted it was found that the supersonic nozzle 
flowfields were typically near chemical equilibrium conditions. This 
is due to the fact that the propellant is typically operating under 
high pressure and high temperature conditions and the nozzle expansion 
ratios are not large. Because vibrational, rotational and translational 
modes of energy storage equilibrate more readily than chemical phenomena 
it is argued that these modes will also be near equil ibri urn conditions. 
This assumption has-been applied throughout this analysis. 



APPENDIX D 


UNSTEADY GAS DYNAMICS EFFECTS 



APPENDIX D NOMENCLATURE 


A = area, 

d = diameter, m 

F = thrust, N 

g = acceleration due to gravity, 9.8 m/s^ 

m = ramjet control volume propellant mass. Kg 

m = propellant mass flow rate, Kg/s 

(mv) = ramjet control volume propellant momentum. Kg m/s 

P = pressure. Pa 

t = time, sec 

V = velocity, m/s 

A = denotes change of a variable 

6 = conical diffuser half-angle, degrees 

o = acceleration, g's 

p = propellant density, Kg/m^ 

T = fluid particle residence time, sec 

A subscript on a variable denotes conditions at the corresponding 
location shown in Fig. D1 . 



APPENDIX D 


UNSTEADY GAS DYNAMICS EFFECTS 

The analyses presented thus far in this thesis have assumed that 
the propellant flow can be considered to be in a quasi-steady state. 

In reality, however, the flowfield variables will be changing with res- 
pect to time as the electrothermal ramjet accelerates down a launch 
tube, and it is possible that errors will be introduced by assuming 
quasi-steady flowfields. It is the purpose of this Appendix to pre- 
sent preliminary calculations which address the effects of unsteadiness. 

In order to analyze unsteadiness effects it is helpful to consider 
the residence time (t) of a typical fluid particle in the ramjet and to 
estimate the flowfield changes that will occur over this residence time. 
Values for t will be calculated from the quasi-steady analyses used to 
generate the results presented earlier in this thesis. Values for fluid 
particle residence times are generated by specifying a geometric con- 
figuration and calculating the time that it takes a fluid particle to 
pass through this configuration when the internal propellant velocities 
relative to the projectile are known. The geometric configuration 
selected (Figure D1 ) has a conical half-angle (6) of 30", an annular 
diffuser throat diameter (d 3 ) of 16 cm, an annular heat addition region 
diameter (d^) of 8 cm, a heat addition region length of 1.0 m, and a total 
projectile length of 1,25 m. Three separate projectile velocity levels 
of 4000 m/s, 7000 m/s and 15,000 m/s will be examined and one value for 
the thrust per unit tube area (F/Aj) will be examined for each velocity. 
For 4000 m/s the thrust per unit tube area is 6.7 x 10’ Pa and for 
7000 m/s and 15000 m/s the values for the thrust per unit tube area are 




rs»a®'»a'S%"«Ta’ f 
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9.5 X 10^ Pa and 1.1 x 10** Pa respectively. The calculated fluid par- 
ticle residence times for 4000 m/s, 7000 m/s and 15000 in/s are 
respectively 


-3 

^4000 ■ 0.769 X 10 sec 
T 7000 = 1 .36 X 10" sec 

^15000 = 0.861 X 10'^ sec. 

The values for the non-dimensional thrust that correspond to these three 
velocities are 

(b-x-) = 6 

4000 


^ Pi ^7000 


= 10 


^'^ 1 ^ 1^5000 


70. 


It is noted that the fluid particle residence times are primarily dom- 
inated by the length of time that a particle spends in the heat addition 
region of a ramjet, and in order for the residence times to be calcu- 
lated with high degrees of accuracy it is necessary that the details of 
the heat addition be known. Also, increasing the amount of heat added 
per unit mass of propellant increases the average propellant velocity 
in the heat addition region. This in turn causes the heat addition 
region residence time to decrease. This explains why there is no 
definite trend for the residence times presented earlier. The differ- 
ences are due primarily to the selection of the heat addition rate per 
unit mass of propellant, and these residence times were selected as being 
representative of what could occur in an actual device operating with 
the geometric configuration specified. 


V 
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During a residence time (t) the ramjet velocity (v) will increase by 
an amount Av. It is instructive to calculate the ratio Av/v, which is 
the ratio of the change in velocity over the time t to the velocity of 
the ramjet at the beginning of the time t. For constant acceleration 
this ratio is given by 

AV og T 

V V 


Using an acceleration level (a) of 30,000 g's, which will be used 
throughout this Appendix, along with the values for the residence times 
just calculated gives the following normalized changes in projectile 
velocity 

=0.057 

'' 4000 


= 0.0t,7 

7000 


'' 15000 " 


These values are reasonable and suggest that the projectile velocity 
will not change significantly during typical fluid particle residence 
times for an acceleration level of 30,000 g's. 

It is also possible that propellant mass and propellant momentum 
will accumulate within the ramjet control volume, and it is instructive 
to consider the control volume flowfield changes over the fluid particle 
residence times. Consider first the possibility of mass accumulation 
for a projectile undergoing a constant acceleration. At any time t the 
influx of mass is given by 
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"i-jp “ AiPiVj . 

The change in the ramjet control volume mass due to the change in the 
incoming mass flow as the projectile accelerates during the residence 
time T can be approximated by 

Am = (m^+^ - m^) T , 

and this can be divided by the mass of propellant in the ramjet (m) to 
determine the ratio Am/m which should be small. 

At any time t the momentum flow entering a ramjet control volume 
is given by 

(fflv) = piAiV^ , 

where (mv) is the ramjet control volume propellant momentum relativf to 
the ramjet, and the amount of momentum that will enter a ramjet control 
volume over a residence time x can be approximated by 

A(«iv) . f"*), ) '■ 

Again the ratio of this momentum change to the propellant momentum 
relative to the ramjet (mv) should be small if unsteady effects are to 
be negligible. Computation of the ratios Am/m and A(mv)/(mv) have been 
performed and the results are shot^n below. For each calculation, it 
has been assumed that the density of the incoming propellant does not 
change over a residence time (though increases in velocity would 
ordinarily imply reductions in density). 
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(M) 


7000 


= 5.57 X 10"^ 


(— ) = 1,32 X 10“^ 

15000 



0.11 


0.45 


0.14 



The above results suggest the errors associated with the changes in 
mass are not large, while the errors associated with the changes in pro- 
pellant control volume momentum are more significant. It is noted, 
however, that proper control of the propellant density immediately up- 
stream of the accelerating ramjet will allow these momentum errors to be 
reduced. The errors that were estimated are felt to be at the 
upper limit of what might actually occur and although a complete analysis 
has not been carried out, intuitively it is felt that the ramjet opera- 
tional conditions could be controlled in a way that would make these 
errors insignificant. For example, it is apparent that the propellant 
density should decrease as ramjet velocity increases so that mass and 
momentum accumulation will be minimized, and this decrease in density is 
also required ior optimum quasi-steady operation (see Figure 12). 

Another area of concern is associated with the fact that the pro- 
pellant in the heat addition region is being accelerated, and it is 
possible that a "hydrostatic" pressure gradient could be produced in 



the heat addition region. This problem was also analyzed for the cases 
presented earlier and results for the ratio aP/P were computed. In this 
ratio aP is the “hydrostatic" pressure difference across the heat addi- 
tion region and P is the average quasi -steady pressure calculated for 
the heat addition region. Computed values of this ratio are 

(^) = 0.064 

4000 

(^) = 0.028 

7000 

(~) = 0.0075 . 

15000 

These values are considered not large enough to introduce significant 
errors into the results presented thus far. 

In conclusion, it is felt that the errors incurred by assuming 
quasi-steady operation, while significant, are not sufficiently large to 
suggest the results presented are grossly in error. Further it is sug- 
gested that these errors can be reduced significantly by proper launch 
trajectory design. It is also noted that these results were obtained 
for a specific geometric configuration and acceleration level. De- 
creasing the various dimensions presented earlier and decreasing the 
acceleration level will decrease the errors associated with unsteady 
gas dynamics effects. 



End of Document 



